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Abstract —Consideration of the nuclear electric quadruple terms in the expression for the fusion 
coulomb barrier suggests that this electrostatic barrier may be substantially modified from that 
calculated under the usual plasma assumption that the nuclei are electric monopoles. This effect is a 
result of the nonspherical potential shape and the spatial quantization of the nuclear spins of the fully 
stripped ions in the presence of a magnetic field. For monopole-quadrupole fuel cycles like p- n B, the 
fusion cross-section may be substantially increased at low energies if the protons are injected at a small 
angle relative to the confining magnetic field. 

1 + INTRODUCTION 

Many nuclei possess nonspherical electric charge distributions* observable as 
electric quadrupole moments (e.g, 2 D, 6 Li, n B etc.). To our knowledge the effects 
of such charge distributions on the physics of thermonuclear plasmas has not been 
studied. We present here some results of simple illustrative calculations for 
several species, at plasma ion temperatures Tj — 40 keV-1 MeV, which suggest 
that such studies may be warranted. 

In the quantum mechanical partial wave model for nuclear interactions, the 
interaction barrier e is taken to be the sum of the nuclear and electrostatic 
potentials at the ‘fusion distance’ d , where d is defined by {de{r)ldr) r=d = 0, An 
analysis of nuclear shape effects on total reaction cross sections at high energy has 
been given by Wong (1973) for s-waves, and of course the role of nuclear 
multipole electric moments in coulomb citation has been studied in detail 
(Alder et al 1956; Bohr* 1976), We discuss here the effect of a nonspherical 
coulomb barrier on the fusion cross sections for Light nuclei, of interest to the 
controlled fusion community, when the collision energy is less than the classical 
fusion barrier height (i.e. the monopole value). In our analysis, we shall define the 
fusion distance d as the distance of closest approach of the two interacting nuclei 
before the nuclear attractive forces come into play. Once the nuclei have 
approached to within this distance* we may assume that fusion will occur. Thus, 
for our purpose, we do not need to consider the nature of nuclear attractive 
potentials* but we need simply to ask what the coulomb potential energy is at the 
distance d. Although a quantum mechanical approach to this problem would be 
somewhat more exact, we use classical physics to illustrate the salient features of 
the problem and to suggest possible areas for future research. 

Our calculations are in three parts: first, we calculate the change in the 
electrostatic potential barriers between selected nuclei at the fusion distance when 
the quadrupole terms are added to the monopole repulsion; secondly, the 
corresponding changes in the fusion ctoss sections are estimated for the most 
desirable (i.e. lowest energy) relative orientations of the fusing nuclei; and finally, 
we show that for quadrupole-monopole cycles (e*g, p- n B) the average fusion 
cross-section is largest when the monopoles are injected at a small angle relative 
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to the confining magnetic field in the containment device. This effect may be 
useful in the start-up of advanced fuel, controlled fusion reactors fuelled by 
neutral injection. 


2. CALCULATIONS OF THE COULOMB BARRIER 
Figure 1 depicts the nuclear model used. Nuclei A and B are separated by the 
fusion distance d AB = a A + a B > where the nuclear effective radius is a = a 0 A 1/3 fm, 
A being the mass number and a 0 — 1.41 fm. When the nuclei are a distance d AB 
apart, the coulomb interaction still predominates (Clayton, 1968 ), At shorter 
distances, strong attractive forces arise. In case (a) the nucleus A has its quad- 
rupole axis perpendicular to the radius vector between the nuclei (r 22 ); ]n case (b) 
this axis is aligned with r 22 . The angles 0, iff thus describe all possible relative 
orientations between the nuclei. In general terms, the ratio of the quadrupole- 
quadrupoie energy € q _^ of the charge configuration at this distance to the 
monopole-monopole energy is simply 


€ q-q _ fl A F fl B y / y 

Z A Z B e 2 A 


( 1 ) 


The monopole-quadrupole energy e m _ q is defined in the same fashion. Figure 2 
shows the calculated variation of this ratio with changing 0, <p, and $ for the 
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Fig. 1.—Geometry for calculation of quadrupole-quadrupole interaction energy for two 
nuclei in 'physical contact’ (separation between centers - + The quadmpole 

moments are Q A = Z A e{a A f ) 2 and Q a = Z H e(a R ') 2 , In case (a) nucleus A has its 0 A 
normal to the line r 22 between centers. In case (b) Q A is aligned with r 22 . 














Fig. 2.—Calculated classical electrostatic energy ratio at the fusion distance as 

a function of orientation angles 8 t <f> and $ for fa) 2 P- 2 D and (b) 2 D- 7 Li. 


2 D- 2 D and 2 D“ 7 Li interactions. (The quadrupoie values have been taken from the 
American Institute of Physics Handbook , 3rd edition, McGraw-Hill, New York, 
1972), Similar calculations for the 2 D- 6 Li case give maximum and minimum 
values of € q _ q /€ m _ m of 1.03 and 0.97. For monopole-quadrupole interactions, the 
modification to the interaction energy at the fusion distance is usually much 
smaller than that of Fig. 2, although some such reactions are substantially 
modified as shown by Table L Of particular interest as an advanced fusion fuel 
cycle is the neutronless p- n B reaction which, although it has a high 'ignition 
temperature 1 , T^SOOkeV or so, offers certain potential advantages over other 
fuel cycles (Dawson, 1976). The nuclear magnetic dipole-dipole interaction (he. 

Table 1.—Classical ranges of the 
ENERGY RATIO AT THE FUSION 

DISTANCE FOR SOME SELECTED 
INTERACTIONS 


Interaction 

Em-m. 


2 D-p 

452 

0,97-1.06 

2 D- 4 T 

378 

0.98-1.04 

2 D- s Hc 

757 

0.98-1.04 

7 Li-p 

1053 

0.80-1.62 

7 Li- 3 T 

914 

0.84-1.43 

**B-p 

1584 

0.83-1.48 












interaction between the spins of the nuclei) energy at fusion distances also 
modifies the total interaction energy, but such effects are small (we find approxi¬ 
mately a 1% modification to e,„., n for 2 D- 2 D due to the magnetic spin interaction; 
note that the spin, magnetic dipole and electric quadrupole axes are colinear). The 
quadrupole deformation of the nuclei will also influence the nuclear potential— 
but this effect is not significant for the present discussion for the reason stated 
earlier. 

It is clear from the simple classical calculations that the coulomb fusion barrier 
is dependent on the relative orientation of the participating nuclei {i.e. fully- 
stripped ions in a fusion plasma) if at least one of the nuclei has an electric 
quadrupole moment. We estimate the corresponding dependencies in the fusion 
cross-sections cr m _ m , and o- q _<, as follows. The effect of the classical potential 
barrier on cr appears in the Gamow factor (Galin el ai, 1974; Glas and MosEt, 
1974), which for convenience we write as 

<r m _ m exp (-0.03 j3e m - m £ -1 ' 2 ) (2) 

where 

and where E is the relative center of mass kinetic energy, A a and A B are the 
respective mass numbers, and the energy units are keV, We estimate cr q „ q and 
a m _ q by simply substituting and e ni _ q for e m _ m in (2). The error involved in so 
doing is not as bad as one might fear, since the form of the interaction potential 
remains substantially the same when the quadrupole terms are included except at 
very close encounter distances. The principal modification to e m _ m is a classical 
lowering or raising of the electrostatic potential barrier height, rather than an 
altered semi-classical barrier penetration probability (i.e. tunneling is weakly 
affected). 

For the two quadrupole-quadrupole interactions studied, we find the range in 
cr q ycr m _ m to be 1.21-0,62 for 2 D- 2 D at E = 40.5 keV, and 1.79-0.002 for 2 D- 7 Li 
at E = 96.6keV. These ranges for the monopole-qnadrupole cases presented 
earlier are shown in Table 2, Although we are mindful of the fact that fusion cross 
sections are not experimentally determined to better than a few percent, the 
indicated 700%-plus enhancement in tr for p- u B at E= 158*4 keV in the most 


Table 2, —Classical estimates of ranges in the fusion 

CROSS SECTION RATIO AT TWO ENERGIES E FOR 

SOME SELECTED INTERACTIONS 


Interaction 

keV 

in 

(£=0.1 e m . m ) 

■ I *ni m 

<E = e m - m ) 

2 D-p 

452 

1.12-0.80 

1.04-0.93 

J D-T 

378 

1.12-0.80 

1.04-0,93 

J D- 3 He 

757 

1.17-0.73 

1.05-0,91 

7 Li-p 

1053 

5.34-0.01 

1.70-0,19 

7 Li- 3 T 

914 

9.30-0.0025 

2.02-0.15 

“B-p 

1584 

7.27-0.0037 

1.87-0.17 






favorable orientation over o- m „ m at the same energy is quite significant. For this 
reaction, er m ^i peaks when 0 = 0 (Fig. la) where particle A is the proton. 

3. FUSION ENERGY BARRIERS IN A MAGNETIZED PLASMA 
Collisions at E :£ e m _ m for the examples above will be principally of three 
kinds: (1) Rutherford scattering without coulomb excitation, since the excited 
states (including rotational bands) in general need collisional energies of a few 
MeV; (2) nuclear potential scattering; and (3) fusion reactions. Although a full 
quantum mechanical treatment of the problem is justified, the classical approach 
used here amply illustrates the important features of the problem and avoids 
unnecessary mathematical obscuration. Only head-on fusion collisions are of 
interest here, so scattering is not discussed. We have chosen the p-"B reaction to 
illustrate how advantage can be taken of the nonspherical e m . q shape; the rest of 
Table 2 is equally amenable to such study. 

A hot n B plasma with mean ion energy (E) = kT; = 150 keV in a magnetic 
field H with ground state nuclear spins f=§, will have ion spatial spin quantiza¬ 
tions m } =|, |, -§ and corresponding energy levels L/ = -g I m I ft|H|/2m p c 

where gj = the nuclear g factor, m p =the proton rest mass, h = 27 rfi = the Planck 
constant, and c = the velocity of light. Frequent collisions, ( E )» U, and a required 
rotational energy change AE R =4.8MeV for the 1 = §-»§ transition ensure that 
the n B nuclear spin states remain unchanged until fusion occurs, and that all the 
spatial orientation states are equally populated. 

Assume that the protons and the n B ions are injected into the containment 
device, and that the boron ions are stripped of all the electrons. Protons are 
injected at angle x relative to H, as shown in Fig. 3. The conical surface is traced 



Fig. 3.—Model geometry for a monopole-quadrupole interaction (e.g. p- n B) in the 
magnetic field H. The angle a ; is quantized with values cos -1 (m,/(I{I+1)) 1 ^) where I is 
the nuclear spin and m, is the magnetic quantum number. 





out by the spin / as it processes through angle «; a, is one of the quantization 
angles, and £ is the angle between the normal to the injected p-beam and the 
quadrupole spin (f = 6 in Fig, 1(a) where particle A is the proton). 

The average value of e m _ m at the fusion distance for all a s = 
cos 1 (m,/(/(7 +1)) 1 ' 2 ) as a function of the injection angle x is given by: 


( q/^m ni) v 


=—i-r 

* ie m-m / = ] Afj J^.(, 


j = tt) 

0) 


(*, 0 df 


( 3 ) 


where A^ = f J (w=7r)-f / (w = 0) I £= tt/2-(x + a r cos m) and n = 2for n B, Thus 
the averaged-over-all angles, normalized {to c m _J interaction energy or fusion 
potential barrier (e) for an ensemble of U B nuclei in a magnetic field is 

2 f^ 2 

— ~ , (4) 
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Numerical integration of (3) gives the curves shown in Fig, 4, The middle 
dashed curve is the arithmatic average of the upper and lower curves, and is the 
representative variation of interaction energy with x since all the spatially 
quantized 11 B spin directions should be equally populated in a high-temperature 
plasma. This calculation also gives {e}^ 1.11, so that the spatially averaged fusion 
barrier for n B is increased by some 11% or so over that expected for the simple 
monopole potential barrier case (i.e. € m _ m for p- n B). 

Figure 4 shows that the most energetically favorable injection angle is ^=0, 
i.e. injection along the magnetic field lines, if advantage is to be taken of the 
quadrupole shape in the electrostatic potential barrier. This advantage is on the 
order of (1.11-0,95/1.1 l)x 100 — 14,4% in barrier height at # = (), diminishes to 
zero at approximately * = 45°, and is a penalty on the order of 11% at # = 90°. 
Note that no appreciable advantage can be gained by polarizing the injected 



Fig. 4— Calculated variation of the classical monopole-quadrupolc interaction energy 
e m^\ } or P“ M E* normalized to € in _ m and averaged over the precession angle w, with 
injection angle x relative to H. The energy (s) is the normalized averaged over all x 

value of 







proton beam, since, ( 1 ) as shown above, the spin-spin (magnetic) interaction 
energy is only - 1 % of € m _ m , and ( 2 ) collisions will rapidly destroy the beam 
polarization as the protons move through the hot plasma. 

4. DISCUSSION 

It is possible to estimate the corresponding changes in the fusion cross section 
tr m _ q for p- n B by using (2). It must be emphasized that in so doing, only an 
approximate value for the cross section is obtained, since we assume that 
(j((€))^<o-{e)> when the results of the previous section are used to estimate 
The cross-section value obtained is an overestimate, but it serves to correctly 
illustrate the dependence of the quadrupole effect on the injection angle. 

For the case of parallel p injection (x = 0) and a relative center of mass 
collision energy of £ = 158keV, a m _ q ^I 0 4 - 10 5 (<r m _ q ) 7 where < 0 —}= ^({€», 
This laTge enhancement is a consequence of the exponential behavior of the 
classical Gamow probability, and should not be taken explicitly. Rather, this is an 
indication of the rapidly changing fusion cross-section in this energy range. This 
means that significantly lower injection energies could be used to initiate the same 
fusion reaction rate obtained at higher energy, if the injection is along H. 

It is not clear whether this quadrupole effect will be generally advantageous 
for steady-state fuel injection, due to plasma physics problems such as ion 
beam-driven instabilities (Mynick et al , 1977; Yamada et al , 1977) and particle 
loss processes (Baldwin, 1977) which occur for injection parallel to the confining 
magnetic field. However, there may still be sufficient energy advantage at small 
injection angles (say 10 °<^< 20 °) to justify further study of this effect as an aid 
in the ignition phase of fusion reactors fuelled by neutral injection. 


5. SOURCES OF UNCERTAINTY 

The separation distance d = a 0 (A + A 3 3 ) at which we evaluate the classical 
electrostatic contribution of the quadrupole moment should be that minimum 
distance inward of which the nuclear attractive forces come into play. Thus any 
uncertainty A a Q in our knowledge of a 0 (he. a 0 = 1.41±Aa 0 fm) enters into the 
discussion of the energy ratios. Similarly, any uncertainty A a r in the quadrupole 
parameters a f (see Fig. 1) affects the energy ratios. Using (1), a little algebra 
shows that the maximum compounded uncertainty in the energy ratios is 


A(e q - q /g m —m) 1 T^ I^AqI . 0 A# | ~| 

€ q-ql € m-m L *0 «' J 

A(€ m - q /e m - m ) . T |Aa 0 | t lAa'H 

m-m L flfi O* J 


( 5 ) 

( 6 ) 


We can now obtain some tentative numerical estimates. For [Aa 0 |/a 0? we use a 
liberal error of 5%. For |A« r |/£*', we use 7 Li as a typical example (Green, 1971} 
and find |Aa'|/a f ^2%. These figures lead to uncertainties in the energy Tatios of 
— ± 14% in (5) and ±12% in ( 6 ). Referring now to Fig. 4 for instance, we see that 
our suggestion of the importance of the quadrupole effect still remains valid: at x = 
0.96 ±0.11 for all states (spatially averaged), and 
1.26 ±0.15 at #-90°. That is, the effective energy barrier should vary by at 






least 3.7%, and possibly by as much as 70%, as the injection angle changes from 
0° to 90°. 

Quantum mechanical treatment of the problem should lead to refined quan¬ 
titative estimates of the various parameters—but inasmuch as we have performed 
our calculations at the distance d where classical physics still remains valid, such a 
treatment cannot alter the main conclusions of this paper—namely that (i) 
quadrupole moments in some cases can substantially modify the fusion barrier, 
and (ii) due to space quantization of the nuclei, such a modification can lead to a 
lower ignition temperature than the one obtained on the basis of the plasma 
physics assumption that all nuclei are electric monopoles. 
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